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Abstract 

We discuss the heavy quark production cross section near partonic 
threshold in hadronic colhsions, including the resummation of leading 
and next-to-leading logarithms arising from soft gluon emission. We 
show how to handle the complications due to the non-universal non- 
leading logarithms. We give analytical results for the qq partonic 
subprocess and numerical results in the DIS scheme for top quark 
production at the Fermilab Tevatron where the qq channel dominates. 
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1 Introduction 



Heavy quark production has been a topic of intense interest in the last few 
years, particularly since the discovery of the top quark at the Fermilab Teva- 
tron. The top quark production cross section will be measured with in- 
creased precision as Tevatron running continues. It is therefore important to 
make precise theoretical predictions for the cross section. Calculations of the 
cross sections for processes such as top quark production are based on the 
factorization theorems of perturbative quantum chromodynamics (pQCD) 
[||]. Factorization separates the perturbative, short-distance hard scattering 
from nonperturbative universal parton distribution functions. The hadronic 
cross section is then given by the convolution of the perturbatively calculable 
partonic cross section with the experimentally determined parton densities. 
Complete calculations of heavy quark production have been carried out up 
to next-to-leading order (NLO) in the strong coupling constant @ and 
the NLO corrections were found to be significant. Calculations to higher or- 
ders are formidable and full results do not exist beyond NLO. However, near 
the final-state production threshold, z = Q"^ / s = 1, where Q is the invari- 
ant mass of the heavy quark pair and s is the square of the center-of-mass 
energy of the partonic collision, there are large logarithms at each order in 
the perturbative expansion originating from soft gluon emission which can 
be resummed to all orders in pQCD. At ?7,th order in perturbative QCD one 
encounters terms as singular as (— a"/n!) [ln^"^^((l — z)^^)/ (1 — z)]+, which, 
when folded with parton distributions, give large and positive corrections [§]. 
These Sudakov logarithms, arising from the incomplete cancelation between 
real and virtual corrections near threshold where the energy of the radiated 
gluon approaches zero, increase the NLO cross section. 

Resummation is a result of factorization 0] and is most easily derived in 
terms of moments of the cross section with respect to the variable r = Q"^/ S, 
where 5* is the center-of-mass energy of the incoming hadrons. In moment 
space the hadronic cross section becomes the product of moments of the 
parton distributions and the partonic cross section. The physical resummed 
cross section is obtained upon inversion of the exponentiated moments back 
to momentum space. 

The leading logarithms for heavy quark production arise from soft gluon 
emission from the quarks and gluons in the incoming hadrons and are there- 
fore universal. They are thus the same as in the Drell-Yan process where 
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Sudakov resummation was applied some time ago ^. More recently, 
leading-log resummed calculations for heavy quark production |0, H, H, |10|, |ll 



have been presented, using different methods [0, ^ |TT| to invert the moment 
calculation to momentum space. The differences between the results are at 
the level of subleading logarithms and are numerically small for top quark 
production at the Tevatron. 

The theoretical framework for the extension of resummation to non- 
universal next-to-leading logarithms (NLL) for heavy quark hadroproduction 
has been presented in Refs. |P, |12|, in moment space, and has recently 
been extended to dijet production and single-particle inclusive cross 
sections [0, including direct photon production |jl5| and heavy quark elec- 
troproduction The formalism has also been used to study dijet rapidity 
gaps Recently a different method for the calculation of the total NLL 
resummed heavy quark cross section has appeared [|TB| which agrees with the 
moment space results in Refs. ||, |T2|, [1^, and has also been applied to direct 
photon production [|Tgf| . 

Beyond leading logarithms, the color exchange in the hard scattering 
must be taken into account. In the next section we will present an anoma- 

which controls 
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lous dimension matrix, explicitly calculated in Refs. 
color-sensitive gluon radiation into the final state. Using the results in 
we have previously calculated the NLL resummed hadronic cross sections for 
top and bottom quark production at fixed center-of-mass scattering angle 
ig. In this paper, we present analytical results for the angle-integrated 



NLL resummed cross section in the qq QQ channel in the DIS scheme. 
The calculation of gg —>■ QQ is similar but more involved and will be dis- 
cussed in detail elsewhere In Section 3, we give numerical results for top 
quark production at the Fermilab Tevatron in the qq —>■ tt channel, which is 
dominant. 



2 Resummation formalism 

The resummation of the cross section in moment space is achieved by refac- 
torizing the partonic cross section into hard components which describe the 
truly short-distance hard-scattering, center-of-mass distributions associated 
with gluons collinear to the incoming partons, and a soft function associated 
with non-coUinear soft gluons. The soft function is a matrix in color space 
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which satisfies a renormahzation group equation whose solution provides a 
matrix evolution equation, in terms of soft anomalous dimension matrices, 
that controls threshold logarithms |]T2|. The relevant anomalous dimension 
matrices were calculated in Refs. [|1^, for heavy quark production through 
light quark annihilation and gluon fusion. At the level of NLL we can diago- 
nalize the matrix equation and calculate the eigenvalues and eigenvectors of 
the anomalous dimension matrix. The resummed cross section then becomes 
a sum of exponentials. 

As mentioned above, we have previously applied these results to top and 
bottom quark production near partonic threshold in hadronic collisions at a 
center-of-mass scattering angle 9 = 90° where the soft anomalous dimension 
matrix is diagonal . Here we consider the general case of the total angle- 
integrated cross section where we must diagonalize the matrix and decompose 
the Born cross section in a particular color tensor basis. The complications 
are, as we will see below, considerable, even for the qq channel. 

We now present an analysis of the process q{pa) + q{Pb) Qipi) + Q(P2)- 
We define a variable S4 = s + ti + Mi in terms of the Mandelstam invariants 
■5 = (Pa + ^1 = {Pa — Vi)"^ ~ "^^j and Ui = {ph — pi)^ — w?, where m is 
the heavy quark mass. The value of S4 depends on the four-momentum k of 



the radiated gluon and is given near threshold by S4 = 2m (1 
The qq resummed partonic cross section is then given by R, 0, 



- 
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2mA;" 



E 



dcosO 



dS4fqq,ij{S4,0) ^ 



dsA 



(2.1) 



where ij is the component in color space of the anomalous dimension matrix, 
and da^ij{s,S4,6)/ds4 are components of the differential of the Born cross 
section defined in the same color basis as the anomalous dimension matrix, 
as explained below. 

The function fqq^ij is given at NLL by the exponential 



fqq,ij = exp[E, 



eXp[Egq + Eqq{Xi, Xj)] 



(2.2) 



where Xij are the eigenvalues of the soft anomalous dimension matrix, Ts- 
The universal component. Egg, is scheme dependent and is known from Drell- 
Yan production P, At NLL in the DIS scheme this universal contribution 
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is given by ^0 



E. 
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1 dto' 



Zn 



(«s(0 

3 Cp /tu'/i^ \ -| 



(2.3) 



2 TT A2 

where Cp = {N^-1)/{2N), Ca = N,K = CA(67/18-7rV6)-5n//9 H, A is 
the QCD scale parameter, and ujq = Si/{2m?) with the number of colors 
and Uf the number of flavors. The NLL color-dependent process-specific 
contribution to the exponent is pO 



Eqq{Xi, Xj 



1 dto' 



A,: 



UJ' 



as 



a. 



A2 



e 



(2.4) 



where i,j = 1,2. The S4 dependence appears only in the argument of the cou- 
pling constant. Note that we have introduced a cutoff Scut in the S4 integra- 
tion in Eq. ( |2.1| ) because, in the exponents, a, diverges when uj'^ fi"^ / A"^ ~ 1, 
corresponding to a minimum S4 of S4,min = 2m? K/ ^. We choose a value of 
the cutoff consistent with the sum of the first few terms in the perturba- 
tive expansion ||^, ||, ||, 0, in the range 30s4^min < s^nt < 40s4^min. Note 
that this corresponds to a cutoff on the soft gluon energy of the order of the 
decay width of the top, giving a natural boundary of the nonperturbative 
region (see also the discussion in |]10|)- The central value in our cutoff range, 
■Scut = 35s4_min, corrcspouds to S4/(2m2) = 0.04 for m = 175 GeV/c^ and 
A = 0.2 GeV. 

The anomalous dimension matrix is calculated at partonic threshold, S4 = 
0, in a color-tensor basis consisting of s-channel singlet and octet exchange. 



Cl — Csinglct — ^ab ^12 ; ^2 — Coctet 

In this basis and in an axial gauge, = 0, the components of Fg are |1T2| 

-C^(L^ + l + 7ri), 



(2.5) 
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r2i = 


r22 = 
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TT I 




TTt 



+ Lg + TTi 



(2.6) 
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where 



1 - 2mVs / 1-/3 



In- 



1 + 



+ Tci 



(2.7) 



with [3 = yl — Am? / s. Note that the matrix Ts is diagonal in this singlet- 
octet basis when /5 = and also when 9 = 90° {ui = ti) for arbitrary (3. It is 
this property that simplifies the calculation of the cross section at ^ = 90°, 
as shown in Ref . pO| . 



The eigenvalues Aj and eigenvectors Cj, with i = 1,2, of the anomalous 
dimension matrix are given, respectively, by 



Ai,: 



Tu + ± (Fn - 122)' + 4ri2r2i 



1/2 



and 



1 



(2.8) 



(2.9) 



for each eigenvalue Aj. Then, if C = (csingiot , Coctet) is the original color basis, 
the diagonal color basis in which we work is C = {c[ , Cj) = CR, where 



ri2 



ri2 



Ai— Fii A2— Fii 
1 1 



R = [ei 62] = 

and the diagonalized anomalous dimension matrix is 

R-^TsR 



(2.10) 



pdiag _ D-ln 



Ai 
A2 



(2.11) 



The diagonal color basis is given explicitly by 

ri2 



c ^ {c[ , 4) = (— 



"Cginglet ~l~ ^Qctet 



r 
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11 A2 — i 11 

or, inversely, C = C'R~^, giving for the octet contribution 

_ (Ai - rn)c; - (A2 - rn)4 

Coctet / \ 

[Al — A2) 



'-'singlet ~l~ ^Qctet 



(2.12) 



(2.13) 



octet 



Since the Born cross section is pure octet exchange, it is proportional to c; 
so that, after squaring Eq. ( |2.13| ), the Born cross section is defined in terms 
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of 14 141^, and . Using the relations c^octet = " l)/4, cl\ng\et = ^'^i 
and Csingiet Coctet = 0, we find 
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n,2l 



Al,2 — Til 



(2.14) 



and 



J J* 



+ 



(Ai-rn)(A2-rn)* ' 4 ' ^^'^^^ 
The resummed partonic cross section can then be written in the diagonal 



basis as 



NLL,res/ 2^ 



2 .1 

d cos 9 



-1 



-2msl/2 



1 rfaJJ^(s,S4,^) 



Scut 



lAi — As 



dsA 



X 



'^rL + |A,-rn|^)e-..- + 



' 4Ar^ 
_ I 



rf2 + |A2-r 
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2 ] Q^qq, 22 



-rl.Re (e^-.-) - 2Re ((Ai - rn)(A2 - rn)*e^ 



(2.16) 



where 



S — S4 



X 



ds4 " " 4s- _ s^)2 _ 4s^2 

(3(s - S4)^ - 8W)(1 + cos^ 6) + 4sm2(l - cos^ 9)] (2.17) 



is the differential of the Born cross section. The explicit expressions for all the 
quantities in Eq. ( |2.16| ) are long but straightforward to derive. For example, 
if we define variables r and by 



at 



(rii-r22)' + 4ri2r2i = ^re 



i<f> 



we can write the eigenvalues as 

1 



Ai,: 



In + r22 ± ^r'/'e'f'/' 



(2.18) 



(2.19) 



7 



Then, we find for the exponentials 



eE,W22) = exp lEqg- — [Rein + ReFaa 

I JuJo LO 



± 



r ' cos ' — 



TT 



where the +(— ) sign for the last term is for Eqg^ii[Eqg^22), and 
Re (e^^^.i^) = exp lEgg - /' ^ [ReFn + ReFas] 1 



X COS 



UJ' 



TC 



(2.20) 



(2.21) 



Note that lm{e^ii'^^) is the same as Iie{e^'''''^^) above but with the cosine 
replaced by the sine of the terms in the last square brackets of Eq. ( 2.21| ). 



3 Numerical results 

In this section we present some numerical results for the exponents and the 
resummed partonic and hadronic top quark production cross sections. 

In Fig. D, the exponents in the resummed cross section, Eg-^ u and -^95,22, 
Eqs. iW^ ) and ( p.20|) , along with the universal contribution. Egg, Eq. 



are given as functions of S4/(2m^) in the DIS scheme with m = 175 GeV/c^, 
y/s = 351 GeV, and A5 = 0.202 GeV to be consistent with the CTEQ 4D 



parton densities Note that the absolute value of Egg^n is given. The 

color-dependent exponents are considerably larger than -Egg,22 at fixed angle. 



as shown in Fig. 1(a) of the first reference of pO|. 

In Fig. H, the partonic top quark cross section, Eq. (|2.16| ), is presented 
as a function oi rj = s/ (4m^) — 1. We show the NLO exact and approximate 
results as well as the NLL resummed result (with Scut = 35s4,min)- The NLO 
approximate cross section is the one-loop expansion of our NLL resummed 
cross section. The lower limit of the t] range of the NLL resummed partonic 
cross section depends on the value of Scut- No cutoff is applied to the NLO 
approximate results. As shown in [0, |1^], the one- loop expansion of the 
resummed cross section agrees analytically near threshold with the NLO 
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approximate results in |2^, and, as is evident from Fig. ^, is an excellent 
approximation of the exact NLO cross section in the threshold region (see 
also Fig. 7 in P5|). We note that the largest contribution to the hadronic 



cross section comes from the region 0.1 < < 1 P, p!3[ 



The NLL resummed hadronic cross section is given by the convolution of 
parton distributions (pi/hi for parton i in hadron with the partonic cross 
section 



NLL,res/ri 2\ J x ( 2\ j. 2\ NLL, res/ n 2 

^^gg.had l-^, m ) = dT \ —(f)q/h,{x,^ )(f)q/hA-,f^ )(Tqq [t b , TTl 

q=u •'TO •'T X X 

(3.1) 



where a^-^^'^'^^lrS, m?) is defined in Eq. ( |2.16| ) and tq = (m+i/m^ + ScntY / S. 

Our numerical results for the tt production cross section at the Fermilab 
Tevatron with \fS = 1.8 TeV are shown in Figs. ^ and ^ as functions of the 
top quark mass. We use the CTEQ 4D DIS parton densities p3|, |2^. Since 
the parton densities are only available at fixed order, the application to a 
resummed cross section introduces some uncertainty. The NLO exact cross 
sections, including the factorization scale dependence, are shown in Fig. ^ 
along with the NLO approximate cross section, calculated with Scut = and 
IJ? = rn?. Again we note the excellent agreement between the NLO exact 
and approximate cross sections. 

The NLL resummed cross section, including the scale dependence, is 
shown in Fig. ^. The scale dependence is significantly reduced relative to 
the NLO cross section. Note that implicit in the change of scale is a change 
of cutoff through the definition of S4 min. To match our results to the exact 
NLO cross section we define the NLL improved cross section 

imp NLL, res NLO,approx , NLO, exact /q q\ 

In Fig. ^ the hadronic improved cross section is shown for /i^ = m? along with 
the variation with Scut- Note that the improved cross section is calculated 
with the same cut applied to the NLO approximate and the NLL resummed 
cross sections. This cut reduces the approximate cross section relative to the 
exact one. The variation of the improved cross section is small over the range 

,min 

< Scut < 40S4 
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4 Conclusions 



We have given explicit results for the resummed heavy quark production cross 
section to next-to-leading logarithmic level. We have presented numerical 
results for the dominant channel, gg, in the DIS scheme, for tt production 
at the Tevatron. At m = 175 GeV/c^ and y/S = 1.8 TeV, the value of the 
improved cross section for qq — * ti with Scut/(2m^) = 0.04 is 5.7 < CTg^^ad < 
6.1 pb for m/2 < fi < 2m, with a central value of 6.0 pb compared to a NLO 
cross section of 4.5 pb at /i = m. At the upgraded Tevatron with = 2 
TeV, the corresponding range is 7.4 < 0"^^^^^^ < 7.9 pb, with a central value 
of 7.8 pb compared to a NLO cross section of 5.9 pb at = m. We find 
that the corrections relative to NLO are larger for larger scales. The gg 
channel is more complicated and a complete analysis will be given elsewhere 
[ PH] . Adding the gg contribution we predict a total cross section of 7 pb 
at = 1.8 TeV, in good agreement with experimental values from CDF, 
a^j = 7.6l};5 pb ||2^, and DO, a^j = 5.5 it 1.8 pb [^. Gluon fusion is more 
important for 6-quark production at HERA-B where threshold resummation 
is also of importance 0. Our formalism can also be naturally extended to 
heavy quark inclusive differential distributions in transverse momentum and 
rapidity M. 
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Figure 1: We show the contributions to the qq DIS exponents, Eqs. (|2.2| ) and 
( |2.20| ), in top quark production for /i = m = 175 GeV/c^ and ^/s = 351 GeV 
as a function of s^/ (2m^). The sohd curve shows the absolute value of Egg u, 
the dashed curve, -^95,22- The dot-dashed curve shows lOEqg, Eq. ( p.3| ). The 
factor of 10 is included to facilitate comparison. 
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0.005 0.010 0.050 0.100 0.500 1.000 

Figure 2: The resummed partonic top quark production cross section, 
Eq. ( p.l6| ), is shown (sohd curve) as a function of ?7 = s/(4m^) — 1 for the qq 
channel in the DIS scheme and fi = m = 175 GeV/c^ with Scut = 35s4_min. 
The NLO exact (dashed) and approximate, with Scut = 0, (dot-dashed) cross 
sections are also given. 
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Figure 3: The NLO exact and approximate hadronic ti production cross 
sections in the qq channel and the DIS scheme are given as functions of top 
quark mass for pp colhsions at the Tevatron energy, = 1.8 TeV. The NLO 
exact cross section is given for /^^ = (sohd curve), 4m^ (lower-dashed) 
and (upper-dashed). The NLO approximate cross section with Scut — 
is shown for //^ = (dot-dashed) . 
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Figure 4: The NLL resummed and improved hadronic ti production cross 
sections in the qq channel and the DIS scheme are given as functions of top 
quark mass for pp colhsions at the Tevatron energy, y/S = 1.8 TeV. The 
NLL resummed cross section, Eq. (p.l|), is shown with Scut = 35s4^min and 
fi'^ = m? (sohd curve), Am? and (both dashed). The NLL improved 

cross section, Eq. ( p.2| ), is given for Scut = 35s4,min (dot-dashed), 30s4_min 
(upper-dotted) and 40s4,inin (lower-dotted). 
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